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Falicov and Kimball proposed a real-axis form for the free energy of the Falicov-Kimball model that 
was modified for the coherent potential approximation by Plischke. Brandt and Mielsch proposed an 
imaginary- axis form for the free energy of the dynamical mean field theory solution of the Falicov- 
Kimball model. It has long been known that these two formulae are numerically equal to each other; 
an explicit derivation showing this equivalence is presented here. 
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The Falicov-Kimball modelB is one of the simplest 
many-body Hamiltonians. It was introduced in 1969 
to describe metal-insulator transitions in a number of 
rare-earth and transition-metal compounds and was 
solved Wrjrtie limit of infinite dimensions by Brand-t and 
MielschoBta. The earlier work of Falicov's groupB was 
modified by PlischkeO for the coherent-potential approx- 
imation to give an explicit formula for the Helmholz free 
energy in terms of integrals over the interacting density 
of states (DOS). Later, Brandt and Mielscha derived an 
exact formula for the Helmholz free energy in terms of 
summations over Matsubara frequencies in the infinite- 
dimensional limit. Numerical evaluation of these two 
forms^fpr the free energy showed that they were indeed 
equalEfo but no explicit derivation of the equivalence has 
appeared. 

We illustrate this equivalence here for the spinless 
version of the Falicov-Kimball model (generalization to 
higher-spin versions .is simple). The spinless Falicov- 
Kimball HamiltonianEJ is 

w = -£ tvcUi + E f J2 flfi + uJ2 4*flfi (i) 



where c\ (ci) creates (destroys) an itinerant electron at 
site i, f} (fi) creates (destroys) a localized electron at 
site i, tij is the Hermitian hopping matrix (which is cho- 
sen to be nonzero only between nearest neighbors), Ef 
is the localized electron site energy, and U is the on- 
site Coulomb interaction between localized and itinerant 
electrons. Chemical potentials p and p,f are employed 
for the intinerant and localized electrons, respectively. 

In the limit where the spatial dimension d becomes 
large, the many-body problem can .be solved exactly 
when the hopping is chosen to scales as t = t* j1\fd. 
In this case, the so-called local approximation becomes 
exact. We sketch the algorithm used to solve the many- 
body problem, in order to establish our notation. 

The local Green's function G(z) can be written as the 



Hilbcrt transform of the noninteracting DOS p(e) 

1 



G(z) = / dep(e) 



z + ji — £(z) — e 



(2) 



with z in the complex plane and S(z) the local self energy. 
Dyson's equation for the local self energy reads 



E(z) =z + n-\{z) -G- 1 {z) 



(3) 



with X(z) the dynamical mean field (which must be de- 
termined self consistently). Solving the atomic problem 
in a time-dependent field yields another equation for the 
local Green's function 



G(z) 



w 



z + (i — X(z) z + (i — X(z) — U 



(4) 



with wo = Zq/Z, wi = Z\j 1 Z (the localized electron den- 
sity), and Z = Z + Z\ (the atomic partition function). 
The symbols Z and Zi can be expressed as infinite prod- 
ucts 



Zq = (1 



A, 



iuj n + n 



(5) 



and 



Zl=e -^/-M/) (1 + e /3(^t/) ) -Q(l 



jJ n + fi — U 

(6) 

where f3 = l/T, and we used the notation A n = X(iu> n ) 
with ioj n = MrT(2n + 1) the fermionic Matsubara fre- 
quency. 

The Brandt-Mielsch form for the Helmholz free energy 

is 

T = -ThxZ -T I dep(e) ln[(iw„ +/x - E n - e)G n ] 
J n 

+ (ifWi + pp c (7) 

with p c the itinerant electron density. Our aim is to 
replace the Matsubara frequency summation of the loga- 
rithmic function by an integral over the real axis. To do 
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FIG. 1: Contours used in various integrals. In panel (a), we 
show the contour C that surrounds all of the fermionic Mat- 
subara frequencies which are indicated by X's. The dotted 
line denotes the real axis. In panel (b), we show the deformed 
contour C' that allows one to replace the integral by one over 
the real axis. 



this we use Eq. (J3J) to write G n = \/{iui n + p, — A„ — E n ) 
and rewrite the sum in Eq. (Q) as 

T ^ M( iuJ n + jU - S n - e)G„] = 



In 1 



In 1 



iuj n + /i — £ 

Xn 



iljJ n + /X - S„ 



(8) 



The function S(z) = l/[z + /x — £(z)] is the irreducible 
part (with respect to the hopping) of the itinerant elec- 
tron Green's function and it possesses the same analytic 
properties as do the Green's functions (a branch cut on 
the real axis with a change in sign of the imaginary part 
above or below the cutt2l). The dynamical mean field 
X(z) also has the same analytic properties. As a result, 
the logarithmic functions in Eq. (ph are analytic func- 



tions above and below the real axis (the only branch 
cut lies on the real axis) and they behave as 1/z for 
\z\ — > oo. This implies that we can express the Matsub- 
ara frequency summation as a contour integral around 
the contour C illustrated in Fig. [j] (a) yielding 

T ki[(iu n + pL - S w -e)G n ] = 



2ni 



dzf(z) 



In 1 



z + (i - £(z) 



In 1 - 



\(z) 



z + p - E(z) 



(9) 



with f[z) = l/[l + exp(/3z)] the Fermi-Dirac distribution. 
The contour C is deformed to C which runs parallel to 
the real axis as shown in Fig. [l] (b). Since there is a 
branch cut on the real axis, the integral over C becomes 
the imaginary part of the integral from — oo to oo 



T = —ThiZ + fifWi+np c 



+ - I doj I dep(e)f(u) 



Imln 1 



-Imln 1 



AH 



lj + /i — E(cj) 



lu + fi — S(w) 
(10) 



Because the sign of the imaginary part of the functions 
that make up the argument of the logarithms is fixed 
above and below the real axis, the value of the imagi- 
nary part of the logarithms is defined to lie in the range 
between —tt and or and n, depending on this sign. 
To satisfy the analytic properties of the logarithms in 
Eq. ( |To| ) , note that the expression in the square brackets 
can be rewritten as 



Im In 



+ fx — £(w) 



uj + /i — E(w) — X(u>) ' 



(11) 



but one must be careful not to shift the imaginary part 
of the logarithm by an integer multiple of 2n, which cor- 
responds to a different sheet of the logarithm. 
Noting that 



/M = -T— ln[l 



•exp(-/3u;)], 



(12) 



allows us to integrate by parts (since the boundary terms 
vanish) and gives 



T 



T 

+ - 

TT 

x Im 



TlnZ 



1 



fXfWl + fip c 

dep(e) ln[l - 



-0U] 



1 - E'(w) - A'(w) 



+ (j, — £(cj) — e lj + /i — X(lj) — E(w) 



with the prime indicating a derivative with respect to lj. 
The integral over e can be performed by using Eq. (j^) 
and the fact that the DOS has unit weight, to yield 



T 



-ThxZ + (XfWi + pp c 
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+ ^ J duj\n[l + e-^}lm[G(uj)X'(uj)}. (14) 

The interacting DOS is defined to be A(ui) = 
— Im[G(w)]/7r. Using this fact, we can add and subtract 
an integral over A(u>) to produce 

+ - [ duj\n[l + e~ l3uJ }lm[G{cj){-l + X'{uj)}}. (15) 



A(w) 



1 - 



UJ + fx — U 



+ (l-wi)ln 



w + fJ, 



Next, we substitue in Eq. (Q) for G(u>) and add 



= ^ J &uln(l + e~ /3 "')lm 
+ T Wl \n(l + e^ l - u) ) 



uj + p-U + i0+ 



(16) 



+ Twiln(l + e /3( ' i - c/) )+T(l- W i)ln(H-e /3M ). (20) 

The contour integral can be evaluated by residues which 
produces a sum over Matsubara frequencies 

T = -T J dujA{uj)\n{l + er f3ul )-T\nZ + p f wi + pp c 

+ ryL^nfi--. — — 

A,, 



and 



1 — Wl 



uj + fi + i0+ 



(17) 



= ^ J duj\n(l + e - /3a, )Im 
+ T(l-Wi)ln(l + e ») 
to Eq. (|l5|). Collecting terms gives 

T = -T J duiA(u)ln(l + er fiu> ) -T\nZ + p fWl + pp, 



+ (l-toi)ln 1- 

V lUJ n + /i, 

+ Tu; 1 ln(l + e /3 ^- c/) )+T(l- Wl )ln(l + e' 3 ^). (21) 

The sum over Matsubara frequencies can replaced by 
terms that involve InZo and IniJi from Eqs. (||) and (|^). 
Collecting terms gives 



T = -T J dujA(uj) ln(l + e-^) 



x lm< 



uj + p - U + i0+ 



(u + H-U)[-l + \'(v)] 



uj + p — U — X(u>) 



1 — Wl 



oj + [i + i0+ 



(w + /i)[-l + A'H] 



uj + p — A(cj) 



+ Tu;iln(l + e /3 ^- c/) )+T(l-u;i)ln(l + e^). (18) 
The terms inside Im{...} can be expressed as a derivative 

T = -T J dujA(uj)\n{l + e- f3u )-T\nZ + p fWl + pp c 
T 



+ Twi In -± + T(l - wi) In + E f w x + pp&Z) 

Using the definitions for wq and wi in terms of the Z's, 
and the relation 



ln(l + e-^) = -/3w/H-/Hln/H 

- [1-/H]ln[l-/H] (23) 

gives us our final result for the Helmholz free energy 
T= \ dujA(uj)f(uj)(uj + (ji) + Eftvi 



+ - / duj\n[l + e~^] 



7T 

x —— Im< Wi In 

doj 



1 - 



*H 



+ (l-wi)ln 



1 - 



uj + p — U 
AH 



UJ + p 



+ Twi\n{\ + e 13 ^-^) + T(l - wi) ln(l + e^). (19) 



Now we integrate by parts and recall Eq. (|12[). Since the 
boundary terms vanish, we are left with an integral over 
the real axis, which can be re-expressed in terms of the 
contour C, and then deformed into an integral over the 
contour C. This gives 

T = -T J dujA(uj)ln(l + e~ 0u ) - ThxZ + fifWi + pp c 



+ Tj duA(uj){f(uj) ln/H + [1 - f(uj)] ln[l - f(uj)]} 
+ T[wilnwi + (l-wi)ln(l-wi)]. (24) 

This is the Falicov-Kimball-Plischke form for the free 
energy which completes the derivation. This form of 
the Helmholz free energy is also correct for the Falicov- 
Kimball model with correlated hopping and it can be 
proved in the same way starting from the expressions of 
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